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Abstract 

We discuss the values of resonance peaks of the cross-section of 
a heavy quark bound state obtained by means of a Green function 
method apphed to a Coulombic model and compare the result to the 
T and J/V' data. 

1 Introduction 

The total cross-section of a heavy quark bound state resonance in an 
e~^e~ annihilation is described by the Breit-Wigner formula 

^(^/^) = M2 _ M)2 + r2/4 

where M is the mass of the resonance, ^/s the centre of mass energy, 
r the total width and Tee the decay width into electrons. 

For a bound state however a better analytical description of its 
total cross-section is given by the imaginary part of the Green func- 
tion of the bound state itself [H El [3] . The basic idea is to consider 
the Schrodinger equation of the bound state and compute its Green 
function 

{-R-E)g{^,y,E) =5{^-y) (2) 
where H is the Hamiltonian of the system 

H=-g + V(x). (3) 
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The imaginary part of the derivative of the Green function given by ([2]) 
taken at the origin is proportional to the cross-section at threshold. 
The finite width of the state is taken into account by the substitution 
E ^ E + iV, E being the energy offset from 2m threshold. So for the 
total cross-section we obtain 

a{E,T) ^ a^+^-^g{<d,Q,E + iV) (4) 

where <Tg+g- accounts for the e+e" part of the process and depends 
on the energy solely by the function 1/s if the energy of the process 
is below the Z threshold. That is, this formula involves only the pho- 
ton propagator and not the Z exchange together with its interference 
effect. The other factors depend on the particular process taken into 
account and are not universal [U [5] . We have therefore factored eq. @ 
in a product of terms, the first one for the e^e~ process and the second 
one for the hadronic part of cross-section. 



2 The Green function method 

The next step is to compute the Green function for a realistic model. 
The Coulombic potential 

V{r) = -^^ (5) 

where r = |x| provides an integrable system which is quite realistic for 
a heavy quark bound state provided the QCD energy scale of is set 
to the inverse of Born radius, rs = 2/(3mag), see [1] and references 
therein. In this case we are able to provide a fully analytical solution 
for the bound state energy levels, namely 

The solution to problem ([2]) in the case of S wave is given in [5] 
with a slight notation change 



TfL 

g{Q,Q,E + iT) = — 



+ log(^) + V(l-i^) + 27-l )] (7) 



where k = — \/m(^ + zr), A = 2asm/3 and the wavenumber u = X/k, 
E = i/s — 2m. The ip is the logarithmic derivative of Euler's Gamma 
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function T{x), 7 ~ 0.57721 is Euler's constant and /i is a soft scale 
that cancels out in the determination of physical observables. 

By inspection of eq. ([7|) we see that the main contribution for ener- 
gies E close to the one of bound states ([6]) is given by the ip function, 
which has simple poles for negative integers —n in the complex plane: 

= —+'ip{n + l) + 

Z + 71 

+ Cik + l)-Cik+l,n+l)]{z + n)'' (8) 

for z — s- — n, where iIj^^^z) is the k-th derivative of the ip{z) function 
with respect to z, the so-called polygamma function; ({z) = J2'k^i 
is the Riemann zeta function, while C{z,q) = + q)~^ is the 

Hurwitz zeta function [6J. 



E 

fc=l L 
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3 Results of the Coulombic model for 

T state 

As previously discussed in sec. [1] the cross-section of a heavy quark 
resonance below Z mass could be written as 

a{E,T) ^ -"^0(0,0, E + iT) (9) 

s 

where in our particular model of Coulombic interaction Q is given 
explicitly by formula d?]). 

With ([9]) we are going to compute the ratio of the 25 and 15 
resonance peaks respectively, which should be independent upon the 
particular bound state chosen for this evaluation. In fact one obtains 
the following expression 

cT{E2s,r2s) ^Mky. ^0(0,0, E2S + iT2s) ..Qs 
a{Eis,ris) 9^(0,0, + iris) ^ ' 

which will be renamed as a{2)/a{l) for sake of brevity. The first term 
of RHS of eq. is close to 1, and from eq. ([6]) the mass of the bound 
state is given by 

M„ = 2m + ^„ = 2m|^l-^^^ (11) 
so that the ratio M^g/Mfg is given hy 1 + a^/S + 0{a^,). 
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Defining a suitable variable depending on a generic width value 
t = T/Ei, with El given by eq. ([6]) for n = 1, we obtain the fairly 
elegant expression 

(t(2) _ 1 (42 7r^ + 425) 
a(T) ~ 8 128 

(147384 C (3) + 1134 TT^ + 11096 tt^ + 36545) 
3072 

+0(t^) (12) 

which would suggest a value close to 1/8 for this ratio. Unfortunately 
the ()12p is a very slowly convergent series as its numeric evaluation 
shows: 

^ = 0.125 + 6.559*2 _ 48,877^4 j_ ^^^6^ /^g^ 

CJ(1) 

that is, the coefficient of 0{t^) term grows more than k\ , suggesting 
to evaluate the result with the full expression of d?]) for the Green 
function. 

We will make use now of this method with T resonances data. 
This bb bound state is the ideal candidate for this method: it is rather 
heavy, mb ~ 5 GeV, so that it should allow us to neglect the confining 
linear term proportional to r in the QCD potential [ZllSlin]. 

For the T(IS') PDG data [TO] give us: 

Mt(IS') = 9.46030 ± 0.00026 GeV , 

r(15) = 54.02 ± 1.25 keV , 
ree(lS) = 1.340 ± 0.018 keV ; (14) 

while for the T{2S) one reads 

Mr (28) = 10.02326 ± 0.00031 GeV , 

T{2S) = 31.98 ± 2.63 keV , 
Tee{2S) = 0.612 ± 0.011 keV . (15) 



Plugging all those data into the Coulombic model ([5]) and eq. ([TO]) 
we compute the ratio of the first two peaks for the T resonance ob- 
taining the value 

4^ = 0.211 ±0.024 . (16) 

a{l) 

This result is very sensitive to the exact total width value, and depends 
much less on both the mass of the state and the exact as coupling of 
the Coulombic model. It differs from first term of the series expan- 
sion (|12p by approximately a factor of 2 proving its slow convergence 
despite the smallness of its parameter t, of the order of 10~^. 
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4 Comparison with data 



We are ready now to compare our results with the T measurements. 
Prom eq. ([T|) the maximal value of the Breit-Wigner resonance cross- 
section is given by the expression 

Using again the experimental data of eqs. (I14|) . (llSp we have the fol- 
lowing value for the measured ratio: 

ai2y 



= 0.685 ± 0.094 , (18) 

' exp 

which is quite sensitive to the exact determination of the widths. It 
is clear that this value doesn't agree with ([160 : some caveats are in 
order here. 

According to [lOlIIl], because of ISR (initial state radiation) and 
Beamstrahlung effect the observed line shape is not simply given by 
eq. ([H) but a convolution of the T width of 0(keV), of ISR and of beam 
energy spread of 0{MeV). The Tee value cannot therefore be directly 
measured, but is calculated from the production cross-section of T 
integrated over the incoming e~^e~ energies, that is 

a(e+e- ^ T)dE (19) 



The integral itself however is again calculated from the Breit-Wigner 
formula ([T]) in a bootstrap fashion, thus leading to a heavily model 
dependent result of the cross-section shape. 

The Coulombic model works well even at lower than T scales [12] , 
there could be however some relativistic non negligible effect in this 
case (recall from sec. [1] that this whole method is non-relativistic) . 
Prom the virial theorem applied to a Coulombic model we obtain the 
relation between the kinetic and potential energies average —{V)/2 = 
(T) that leads to a speed estimate of the component quark inside the 
meson 

{v') = lal (20) 

This brings a 7 relativistic correction of about 8% at T scale that 
could change the computed ratio. 



5 Further estimates: the J/ip case 

To shed some more light on this discrepancy we have described it could 
be useful to compare our results to the one obtained from J/-^ data. 



5 



This cc bound state is less ideal than the former T for our purpose, 
as it is lighter (nic ~ 1.5 GeV), thus linear confining terms for the 
potential as well as relativistic corrections could have larger effects. 
From PDG data for J/^lJ [T^ we read: 

Mj/^(15) = 3096.916 ± 0.011 MeV , 
r{lS) = 93.2 ± 2.1 keV, 
FeeilS) = 5.55 ± 0.14 keV ; (21) 

while for the ip{2S) one reads 

^^(25) (25") = 3686.09 ± 0.04 MeV , 
T{2S) = 317 ± 9 keV , 
Tee{2S) = 2.38 ± 0.04 keV . (22) 

Proceeding like described in sec. Hlwe obtain for the J/ij: peak mea- 
sured values 

= 0.0890 ± 0.0083 . (23) 

The Green function approach method detailed in sections [2] and [3] gives 
us 

44 = 0-0368 ± 0.0019 . (24) 
cr(l) 

A comparison of eq. ()24p with eq. ()23p shows that the result for the 
J /tp case is slightly worse than the one for T as expected. In fact 
for the former case the two central values differ for more than 27 gq, 
while in the latter the difference is less pronounced, above 19(Jg (here 
ac refers to the error given by the Green function procedure). An 
estimate of the relativistic 7 correction for J /'ip done in the same 
fashion of eq. (f20l) gives us a result of about 38%. 

This 7 correction together with the absence of a linear term in the 
potential could account for the larger difference seen at a lower energy 
scale. 

We must stress again that this method of the Green function for 
a Coulombic potential of a Schrodinger equation furnishes us with 
an exact analytical solution given by eq. ([7]) only in the case of this 
particular interaction, and only for a non-relativistic system. The 
addition of further correction terms to the potential of eq. ^ like a 
linear confining term r or a relativistic correction of order would 
spoil the integrability of this problem and thus the possibility of full 
control over solutions. For instance, a Hamiltonian with a funnel 
potential like V{r) = — 4as/ (3r)+ar or a Hamiltonian with relativistic 
correction H = /{2m) — p^/(8m^) + V{r) are not exactly solvable, 



6 



and even a perturbative approach is unfit in our case as corrections to 
the original H of eqs. ([3]) and ([5]) are not small. Those systems would 
call for a purely numerical search of results which is beyond the scope 
of this paper. 

6 Conclusions 

We have compared the ratio of the first two cross-section resonance 
T peaks computed from a QCD model and the one obtained from the 
experiments, given respectively by 0.211 it 0.024 and 0.685 it 0.094. 
Another comparison is done for the J/V' case, giving a calculated 
value of 0.0890 ± 0.0083 and an experimental one of 0.0368 ± 0.0019 
for peaks ratio. Albeit the two results do not agree with each other, 
it is necessary to consider that the theoretical model could need some 
relativistic corrections. On the other hand, the experiments do not 
give a direct measurement of the peaks, but rather depend on the 
model used to evaluate the cross-section. 

Therefore it should be possible to compare again the two results 
using more refined methods from the theoretical model side and from 
the measure technique as well. 
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